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Abstract. The solitary waves of massive (l + l)-dimensional nonlinear S^-sigma models are 
unveiled. It is shown that the solitary waves in these systems are in one-to-one correspon- 
dence with the separatrix trajectories in the repulsive TV-dimensional Neumann mechanical 
problem. There are topological (heteroclinic trajectories) and non-topological (homoclinic 
trajectories) kinks. The stability of some embedded sine-Gordon kinks is discussed by means 
of the direct estimation of the spectra of the second-order fluctuation operators around them, 
whereas the instability of other topological and non-topological kinks is established applying 
the Morse index theorem. 
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1 Introduction 

The existence and the study of solitary waves in (1 + l)-dimensional relativistic field theories 
is of the greatest interest in several scientific domains. We specifically mention high energy 
physics, applied mathematics, condensed matter physics, cosmology, fluid dynamics, semicon- 
ductor physics, etcetera. Solitary waves have been discovered recently in massive nonlinear 
sigma models with quite striking properties, see e.g. [1, 2, 3] and references quoted therein. 
The nonlinear S^-sigma models are important in describing the low energy pion dynamics in 
nuclear physics [4] and/or featuring the long wavelength behavior of 0(N + l)-spin chains in 
ferromagnetic materials [5]. 

In this work, we present a mathematically detailed analysis of the solitary wave manifold 
arising in a massive S^-sigma model that extends and generalize previous work performed 
in [1] and [2] on the N = 2 case. In these papers the masses of the N = 2 meson branches 
are chosen to be different by breaking the O(N) degeneracy of the ground state to a Z2 sub- 
group. Remarkably, we recognized in [1] and [2] that the static field equations of this non- 
isotropic massive nonlinear §> 2 -sigma model are the massive version of the static Landau-Lifshitz 
equations governing the high spin and long wavelength limit of ID ferromagnetic materials. This 
perspective allows to interpret the topological kinks of the model as Bloch and Ising walls that 
form interfaces between ferromagnetic domains, see [6, 7, 8]. 

*This paper is a contribution to the Proceedings of the Eighth International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 21-27, 2009, Kyiv, Ukraine). The full collection is available at 
http:/ /www. emis.de/journals/SIGMA/symmetry2009. html 



2 



A. Alonso Izquierdo, M.A. Gonzalez Leon and M. de la Torre Mayado 



We shall proceed into a two-step mathematical analysis: First, we shall identify the collection 
of N sine-Gordon models embedded in our nonlinear S^-sigma models. The sine-Gordon kinks 
(and the sine-Gordon multi-solitons) are automatically solitary waves (and multi-solitons) in 
the nonlinear sigma model. Second, we shall identify the static field-equations as the Newton 
equations of the repulsive Neumann problem in the S^-sphere. The mechanical analogy in 
the search for solitary waves in (1 + l)-dimensional relativistic field theories, [9], is the re- 
reading of the static field equations as the motion equations of a mechanical system where the 
field theoretical spatial coordinate is re-interpreted as the "mechanical time" and the potential 
energy density is minus the "mechanical potential". This strategy is particularly efficient in 
field theories with a single scalar field because in such a case the mechanical analogue problem 
has only one degree of freedom and one can always integrate the motion equations. 

The mechanical analogue system of field theoretical models with N scalar fields has N degrees 
of freedom and the system of N motion ODE's is rarely integrable. One may analytically find 
the solitary waves in field theories such that their associated mechanical system is completely 
integrable. Many solitary waves have been found in two scalar field models related to the 
Gamier integrable system of two degrees of freedom and other Liouville models, see [10]. The 
extremely rich structure of such varieties of kinks has been exhaustively described in [11]. The 
procedure can be also extended to systems with N = 3 (and higher) degrees of freedom that are 
Hamilton- Jacobi separable using Jacobi elliptic coordinates, see [12]. The mechanical analogue 
system of the non-isotropic massive nonlinear S^-system is the Neumann system, a particle 
constrained to move in a S^-sphere subjected to elastic repulsive non-isotropic forces from the 
origin. It happens that the Neumann system [13, 14, 15] is not only Arnold-Liouville completely 
integrable - there are N constants of motion in involution - but Hamilton-Jacobi separable using 
sphero-conical coordinates. Because the elastic forces are repulsive the North and South poles 
are unstable equilibrium points where heteroclinic and homoclinic trajectories start and end 
providing the solitary waves of the nonlinear sigma model. 

The analysis of the stability of solitary waves against small fluctuations is tantamount to the 
study of the linearized wave equations around a particular traveling wave. We shall distinguish 
between stable and unstable sine-Gordon kinks by means of the determination of the spectra of 
small fluctuations around these embedded kinks. The instability of any other (not sine-Gordon) 
solitary wave in the nonlinear § 3 -sigma model will be shown by computing the Jacobi fields 
between conjugate points along the trajectory and using the Morse index theorem within the 
framework explained in [16]. 

The structure of the article is as follows: In Section 2 we present the fundamental aspects of 
the model. Section 3 is devoted to the study of the set of sine-Gordon models that are embedded 
in the nonlinear §^-sigma model. In Section 4, using sphero-conical coordinates, we show the 
separability of the Hamilton-Jacobi equation that reduces to N uncoupled ODE. The Hamilton 
characteristic function is subsequently found as the sum of N functions, each one depending only 
on one coordinate, and the separation constants are chosen in such a way that the mechanical 
energy is zero. Then, the motion equations are solved by quadratures. The properties and 
structure of the explicit analytic solutions are discussed respectively in Sections 5 and 6 for 
the N = 2 and N = 3 models. Finally, in Section 7 the stability analysis of the solitary wave 
solutions is performed. 

2 The massive nonlinear S^-sigma model 

Let us consider N + 1 real scalar fields arranged in a vector-field: 
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where (y ,?/ 1 ) stand for the standard coordinates in the Minkowski space IR 1,1 , and the fields 
are constrained to be in the S^-sphere, i.e. 

<S 6 Maps(R 1 ' 1 ,S Af ) , = (% + ... + (j) 2 N+1 = R 2 . (2.1) 

The dynamics of the model is governed by the action functional: 



S[$] = J dy*dy l C{d^^) = f dy^dy 1 • - , 



where V(3>) is the quadratic polynomial function in the fields: 

V($) = \ {al4>l + a 2 2 4 + ■■• + a 2 N+1 4> 2 N+1 - a 2 N+l R 2 ) . (2.2) 

We follow standard Minkowski space conventions: y^y^ = g^y^yu, g^ v = diag(l, —1), d^d^ = 
g^d^du = d 2 — d 2 . Working in the natural system of units, K = c = 1, the dimensions of the 
fields and parameters are: [0 a ] = [R] = 1, [a a ] = M, a = 1, . . . , N + 1. 
The field equations of the system are: 



d * am a \ " T^T = - di<t> a = -a z J a + \<f> a , a = l,...,N + l, (2.3) 

^\d(a^ a )J d(pa 

A being the Lagrange multiplier associated to the constraint (2.1). A can be expressed in terms 
of the fields and their first derivatives: 




- do® ■ + &!$ ■ di$ j 

by successive differentiations of the constraint equation (2.1). There exists N + 1 static and 
homogeneous solutions of (2.3) on §> N : 

a = ±R, 06 = 0, Vb^a, X = a 2 a , a = l,...,N + l. (2.4) 

We shall address the maximally anisotropic case: a a ^ a^, Vo / 6, and, without loss of 
generality, we label the parameters in decreasing order: a 2 > a 2 > ■ ■ ■ > ot 2 N+l > 0, in such 
a way that V(®), on the S^-sphere (2.1), is semi-definite positive. 

According to Rajaraman [9] solitary waves (kinks) are non-singular solutions of the field 
equations (2.3) of finite energy such that their energy density has a space-time dependence of 
the form: e{y^,y 1 ) = e(y 1 — vy°), where v is some velocity vector. The energy functional is: 

E[$] = J dy 1 Qflb* • do® + \di$ ■ di$ + V[$) \ = J dy l £{y^y l ) 

and the integrand e(y°,y 1 ) is the energy density. Lorentz invariance of the model implies that 
it suffices to know the y°-independent solutions ^(y 1 ) in order to obtain the traveling waves of 
the model: ^(y ,?/ 1 ) = ^(y 1 — vy°). For static configurations the energy functional reduces to: 

and the PDE system (2.3) becomes the following system of iV + 1 ordinary differential equations: 
' / '"" a 2 J a -\<p a . (2.6) 



1\2 
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The finite energy requirement selects between the static and homogeneous solutions (2.4) 
only two, the set Ai of zeroes (and absolute minima) of V(3>) on S^: 

M = {v + = {o,...,o,R), v- = {o,...,o,-R)}, 

i.e. the "North" and the "South" poles of the S^-sphere. Moreover, this requirement forces 
in (2.5) the asymptotic conditions: 

lim ^-r=0, lim (2.7) 

y 1 -+±oo dy 1 y l ->±oo 

Consequently, the configuration space, the space of finite energy configurations: 

C = {Maps(R,S n )/£ < +00} , 

is the union of four disconnected sectors: C = Cnn U ^ss U ^ns U ^sn, where the different sectors 
are labeled by the element of M. reached by each configuration at y 1 — > —00 and y 1 — > 00. 
Cns and Csn solitary waves will be termed as topological kinks, whereas non-topological kinks 
belong to Cnn or Css- 

3 Embedded sine-Gordon models 

Before of embarking ourselves in the general solution of the ODE system (2.6), we explore 
Rajaraman's trial-orbit method [9]. We search for solitary wave solutions living on special 
meridians of the S^-sphere such that the associated trajectories connect the North and South 
poles, to comply with the asymptotic conditions of finite energy (2.7). In [1] and [2] this strategy 
was tested in the A = 2 case. It was easily proved that the nonlinear sigma model reduces to 
a sine-Gordon model only over the meridians in the intersection of the coordinate planes in M 3 
with the S 2 -sphere, that contain the North and South poles. Consequently, there are sine-Gordon 
kinks (and multi-solitons) embedded in the nonlinear sigma model living in these meridians. The 
generalization to any A is immediate: the A meridians m a , a = 1, . . . , N, determined by the 
intersection of the <f> a — 4>n+i plane with the S^-sphere, all of them containing the points N 
and S, 

m a ^<p b = 0, Vb^a, 6=1,..., A, 4> 2 a + <j? N+1 = R 2 a = l,...,A, 
are good trial orbits on which the field equations (2.3) reduce to the system of two PDE's: 

dl4> a - df4> a = -al4> a + \(/) a , (3.1) 

do^AM-i ~ df&N+i = -alf+i^N+l + (3.2) 

Use of "polar" coordinates in the (j) a — 4>n+i plane: 

4>n+i = RcosO, (fi a = ±Rs'm0, 9 £ [0, ir] 

reveals that the two equations (3.1), (3.2) are tantamount to the sine-Gordon equation in the 
half-meridians sign(</> a ) = ±1: 

d^9-dfe = -l(a 2 a -a 2 N+1 )sm2e. (3.3) 

For the sake of simplicity it is convenient to define non dimensional parameters and coordinates 
in the Minkowski space Mr> : 

°a = ~2 2 ' l = 0- 1 >CJ 2 >--->a A r>U = <T N+1 , 

a(-a z N+1 

af-a 2 N+1 y°, x = J a\ - a 2 N+l y l . 
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Thus, the N sine-Gordon equations (3.3) become: 

2 

□0=($-S£)0 = -^sin 20, a = l,...,N. (3.4) 

The kink/antikink solutions of the sine-Gordon equations (3.4) are the traveling waves 

6^ a (t,x) = 2aTctane ±aMx - vt) - X0 \ 1 = (l - v 2 )^ 

that look at rest in their center of mass system v = 0: 

6% a (x) = 2arctane ±fTa(a; - :Co) . 

The kinks 9^ a belong to the topological sector Cns of the configuration space, whereas the anti- 
kink Oft lives in Csn- Henceforth, we shall term as topological kinks to these solitary waves. 
In the original coordinates in field space, (f> a and cpN+i, & K correspond to two kinks and their 
two anti-kinks, depending of the choice of half-meridian. We will keep the K a notation for the 
4> a > half-meridian and introduce K* for the cj) a < choice: 

R 

<Pa a ( x ) = r 7 n > <I>n+1 ( x ) = ±-Rtanhcr a (x - x ), 

coshcr a (x — Xq) ^ 

LK* _ R ,K* 



4> a a ( x ) = r / s, 4> N a +1 {x) = ±Rt&nha a {x - x ), 

coshcr a (2; — xq) ^ 

<pf a (x) = (pf a (x) = 0, V6^o,JV+l. 



The energy functional, for static configurations, is written, in non-dimensional coordinates as: 

*® = * M \t. 4 +\ e - # (* - e ) ■ m 



where v = yaf — c^+i- The embedded sine-Gordon kink energies are: 
Ex a = Ek* = 2uR 2 a a , Ek x > Ek 2 > • • • > Ek n - 



4 The moduli space of generic solitary waves 

The system of ordinary differential equations (2.6) is the system of motion equations of the 
repulsive Neumann problem with N degrees of freedom if one interprets the coordinate y 1 as 
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the "mechanical time" and the field components <j) a as the coordinates setting the "particle 
position", the so-called mechanical analogy [9]. It is a well known fact that the Neumann 
system is Hamilton-Jacobi separable using sphero-conical coordinates [13, 14, 15]. 

Let us introduce sphero-conical coordinates (u\,U2 . . . ,un) in with respect to the set of 
constants a 2 = 1 — a 2 : 

o\ = < ui < of < u 2 < ■ ■ ■ < u N < a% +1 = 1. (4.1) 
The Cartesian coordinates in 1 JV+1 are given in terms of the sphero-conical coordinates by: 

_ o •> Af Af+1 



' t0 W§!v U(z) = H(z-u c ), B(z)=H(z-a 2 b ). (4.2) 

c=l 6=1 

The algebraic equation characterizing the S^-sphere in R N+1 is simply uq = R 2 in sphero-conical 
coordinates. 

The change of coordinates (4.2) provides a 2 N+1 to 1 map from the S N sphere to the interior of 
the hyper-parallelepiped (in the (ux, . . . ,un) space) Pjv characterized by the inequalities (4.1). 
The map from each "2^ +1 -tant" of the S^-sphere to the interior of Pat is one-to-one. One 
"2 Ar+1 -tant" is the piece of that lives in each of the 2 N+1 parts in which M^" 1 " 1 is divided 
by the coordinate A-hyperplanes. The intersections of the S^-sphere with the coordinate N- 
hyperplanes, N + 1 S Ar_1 -spheres, are 2 N to 1 mapped into the boundary of Pat, a set of 
(N — l)-hyperplanes in the [u\, . . . , «jv) space. Simili modo S^^-spheres, intersection of §> N 
with coordinate (N — k + l)-hyperplanes of K +1 , are 2 N ~ k+1 to 1 mapped into (N — k)- 
hyperplanes of the boundary of Pat. Finally, the North and South poles of S N are mapped into 
a single point in the boundary of Pat: v^ 1 = {u\ = 0, i±2 = of, . . . , un = cr 2 N )- 

The Euclidean metric in ds 2 = d<p\ + • • • + d(f) 2 ^ +l becomes in Pat: 

,2i f , 2 -R 2 u'M fA o, 

ds Is" = }-^ 9aadU a> 9aa = ~4~ B{u ) ' ( ' 

Calculations using sphero-conical coordinates remarkably simplify recalling the Jacobi lemma, 

see e.g. [14]. Let be a collection of p distinct constants (ax, a,2, ■ ■ ■ a p ), ai ^ aj, V i ^ j . Let us 

v 

define the function f(z) = FJ (z — aj), p > 2. The following identities hold: 

i=i 

Now, use of (4.4) allows us easily determine the V(&) function (2.2) in sphero-conical coordinates: 

V(Ux, ...,«„)= 1 - 2^ a a + 2^ Ub } ■ 

\ a=l 6=1 / 

The energy functional for static configurations (3.5) in this system of coordinates can be written 
as follows: 



E\ux, • • • , un] = v / dx < 



Af 

a=l v ' a=l 
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The p — 2 identities on the left in (4.4) allows us to write the potential energy density in different 
ways. We choose an expression (4.5) in the form of a sum of products in the factors (u a — of) 
that guarantees a direct verification of the asymptotic conditions (2.7). 
Alternatively, the energy can be written in the Bogomolnyi form [17]: 

i f dx ( ( du a nn dW\ 2 \ f , dW du a s 

EN ,..., a „] = ,y T y>„(— J +"/^E^^r. («) 

^.a=l x / ) a=l 

where g aa denotes the components of the inverse tensor metric, and W(u\, . . . , un) is a solution 
of the partial differential equation: 

N 

Erg ■ (4-7) 



2 h U'(u a ) "2i-« { au J ■ 



a=l a=l 

Note that (4.7) is the "time" -independent Hamilton-Jacobi equation for zero mechanical energy 
(i\ = 0) of the "repulsive" Neumann problem. Therefore, W is the zero-energy time-independent 
characteristic Hamilton function. 

The ansatz of separation of variables W(u\, . . . , un) = W\(u\) + • • • + Wn(un) reduces 
the PDE (4.7) to a set of N uncoupled ordinary differential equations depending on iV sep- 
aration constants i a . The separation process is not trivial (see [14]) and requires the use of 
equations (4.4). Let write explicitly (4.7) as: 

~ 2B[aa) - {W'M? ~ ~tt TTK - *l) ) = h = 0. (4.8) 





R 2 

The identities of the Jacobi lemma (4.4), for k < N — 2, allow to express (4.8) as a system of N 
equations: 

^|^«K)) 2 - y II K - *D = + nu N a -* + ■■■ + i N . lUa + lN (4.9) 

b=l 

depending oniV-1 arbitrary separation constants 12, ■ ■ ■ , in- 

Asymptotic conditions (2.7) imply, in sphero-conical coordinates, the limits: lim u a = a\. 

x— >-±oo 

This requirement fixes, in (4.9), all the separation constants i a to be zero, that correspond to 
the separatrix trajectories between bounded and unbounded motion in the associated mecha- 
nical system. Thus, integration of (4.9) gives the W function, i.e. the time-independent, zero 
mechanical energy, characteristic Hamilton function of the repulsive Neumann problem: 

N 

W( Ul ,...,u N ) = R 2 Y,(-l) Sa VT^r a , 8 a = 0,1. (4.10) 

a=l 

It is interesting to remark that general integration of the Hamilton-Jacobi equation for the 
Neumann problem involves hyper-elliptic theta- functions [15]. The asymptotic conditions (2.7) 
guaranteeing finite field theoretical energy to the solitary waves - and finite mechanical action 
to the associated trajectories - force theta functions in the boundary of the Riemann surface 
where degenerate to very simple irrational functions. 
The ordinary differential equations: 

du\ _ u dW duM _ nn 

dW 

dx du\ ' dx duN 



8 



A. Alonso Izquierdo, M.A. Gonzalez Leon and M. de la Torre Mayado 




Figure 2. K\ (red) and K2 (blue) kink orbits over S 2 (left). Energy densities for K\jK\ (red) and 
K2/K2 (blue) kinks for a concrete value of < a < 1 (right). 



can alternatively be seen as first-order field equations for which the field theoretical energy 
is of the simple form: E[u±, u%, . . . , itjv] = v J dW, see (4.6), or, as the motion equations of 
the mechanical system restricted to the sub-space of the phase space such that i a = 0, Va = 
1, 2, . . . , N. In differential form they look: 

dx {-l) 5a du a 

- n , a = l,...,N. (4.12) 



b=i 

Although it is possible to integrate equations (4.12) for an arbitrary value of N (see [18] for 
instance), we shall devote the rest of the paper to analyze the explicit solutions and the structure 
of the moduli space of solitary waves in the lower N = 2 and N = 3 cases. 

5 The moduli space of § 2 -solitary waves 
5.1 Embedded sine-Gordon kinks on § 2 

There are two sine-Gordon models embedded in the = and <p2 = meridians, see Fig. 2 
(left). The non-dimensional parameters, for N = 2, are: a\ = 1, a-i = a and 03 = 0. The 
Ki/Kl and K2/K2 sine-Gordon kinks/antikinks profiles read: 

<p? 1 (x) = — — r, <p^(x) = 0, <f)f 1 (x) = ±Rtanh(x-x ), 



cosh(x — xq) 



<P? 2 (x) = 0, <p2 2 (x) = — - — ^ r, 0f 2 (x) = ± J Rtanha(x - x ), (5.1) 

COSh(j(x — Xq) 

(fif 1 (x) = — r, (pf 1 (x) = 0, (ft^ 1 (x) = ±Rt&nh(x — xq), 



COsh(x — Xq) 



-R 



</,; i 2( x ) = 0, (t>2 2 {x) = — : — -P r, <pT 2 (x) = ±Rt<m\iG(x-x ). (5.2) 

COShfj(x — Xq) 

The Ki, K^, K2 and kinks belong to CnSj see Fig. 1; the (^-component is kink shaped for 
the four kinks, whereas the 4>\- and </>2-components are bell shaped respectively for K\ and K2 
but anti-bell shaped for and K%. The antikinks live in Csn and have similar profiles. These 
solitary waves are topological kinks, and they correspond to heteroclinic trajectories in the 
analogous mechanical system. Their energies are: 

E Kl = E K * = 2uR 2 , E K2 = E K * = 2vR 2 a. 

The energy of the K\jK\ kinks is greater than the energy of the K2/K2 kinks, see Fig. 2 (right). 
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5.2 Generic solitary waves on §> 2 

We shall denote (Ai,A 2 ) the sphero-conical coordinates in the 8 2 -sphere, Ao = R 2 , to avoid 
confusion with the N = 3 case. In terms of these coordinates the fields read: 

<Ai = §A!A 2 , ( pl = J^(a 2 -X 1 ){X 2 -a 2 ), 0§ = ^(1 - Ai)(l - A 2 ). (5.3) 

We recall the definition: a 2 = 1 — a 2 , and the (Ai, A 2 ) ranges: 

< Ai < o 2 < A 2 < 1. (5.4) 

The change of coordinates (5.3) maps the § 2 -sphere into the interior of the rectangle P 2 defined 
by the inequalities (5.4) in the (Ai, A 2 )-plane. In fact, the map is eight to one, due to the squares 
that appears in (5.3). Each octant of the S 2 -sphere is mapped one to one into the interior of P 2 , 
whereas the meridians of S 2 contained in the coordinate planes are mapped into the boundary 
of P 2 . The set of zeroes M. of V{(j>\, 2 , 4>z) in S 2 , the North and South poles, becomes a unique 
point in P 2 : = (4>i, 4>2, 03 ) = (0, 0, ±R) =4> = (Ai, A 2 ) = (0, a 2 ). The quadratic polynomial 
V(0i,02)03) becomes linear in sphero-conical coordinates: 

y(A 1 ,A 2 ) = ^(A 1 + A 2 - C T 2 ) 
and the energy functional for static configurations in terms of (Ai, A 2 ) is: 

*m \ i (a J 1 ^ f dX A 2 , R2 f ^i-v 2 ) , A 2 (A 2 -a 2 ) ^\ 

E[\uX 2 ] = vJd X {-J2 9 ^— j Ai _^ + ^_ Ai jj. (5.5) 

Here, the components of the metric tensor (4.3) read: 

)2/\. \_\ E>2 



R 2 (X 1 -X 2 ) -R 2 (X 2 -X 1 ) 



■ yjJ 4A 1 (a 2 -A 1 )(l-Ai)' 522 4A 2 (a 2 -A 2 )(l-A 2 )- 
The knowledge of the solution (4.10) for N = 2: 

W(\ u A 2 ) = R 2 ( (-1) 51 v/l-Ai + (-1) 52 v 7 ! - A 2 ) , «5 2 = 0, 1, 



suggests to write the energy functional (5.5) in the Bogomolnyi form (4.6) leading to the system 
of first-order ODE's (4.11): 

dA i = ( 2A i( A i ~ ^ 2 )v / l^AT dX 2 = S2 2A 2 (A 2 - ^ 2 )yT^A2" 
<ix (Ai — A 2 ) dx (A 2 — Ai) 

These equations can be separated in a easy way: 

= d^i , rfA 2 f56 s 

(-l) 5 i2Ai(Ai - a 2 )VT^AT (-1)*»2A 2 (A 2 - a^VT^' ' 

d2; = dX ^ + rfA2 f 5 7 ) 

(-l) 5l 2(Ai - a 2 )vT^AT ^ (-1)^2(A 2 - a 2 ) v T^A 2 - 1 ' ' 

The Hamilton-Jacobi procedure precisely prescribes the equation (5.6) as the rule satisfied by 
the orbits of zero mechanical energy whereas (5.7) sets the mechanical time schedule of these 
separatrix trajectories between bounded and unbounded motion. 
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Instead of attempting a direct solution of the system (5.6), (5.7) we introduce new variables 
sx = (-l) 5l Vl - AT, s 2 = (-1) 5 V1 - A 2 . In terms of s\ and s 2 the ODE's system (5.6), (5.7), 
after decomposition in simple fractions, becomes: 

E : ^ Sa 9 = dx, 9 = dx. (5.8) 

a=l " a=l " 

Integrating equations (5.8) in terms of the inverse of hyperbolic tangents/cotangents, and using 
the addition formulas for these functions we find the following general solution depending on 
two real integration constants 71 and 72: 

2 ^ 

= h, —, r=*2j ti = tanh(x + 71), t 2 = tanhcr(x + 72). (5.9) 

I + S1S2 <7(Sl+S 2 ) 

To solve the system of equations (5.9) separately in si and S2 we introduce the new variables 
A = s± + S2 and £> = siS2- The subsequent linear system in A and -B and its solution are: 

A- h B = t u *t 2 A-B = o\ A=^^, B = ^f. 

1 — Ol\l2 1 — (Jt\ti 

Therefore, s\ and S2 are the roots of the quadratic equation 

s 2 - As + B = 0, si = , s 2 = - • 

Because Ai = 1 — sf , A2 = 1 — s 2 , the come back to Cartesian coordinates (5.3) provides the 
explicit analytical formulas for the two-parametric family of solitary waves: 

„ sech(x + 7i) 

4> x {x) = (-l) ei Ra- 1 / 1; ei = 0, 1, 

1 — <7tanh(a; + 71) tanncr(x + 72) 

a ( \ ( ^e 2D - tanh(x + 7i)sechcr(x + 7 2 ) 

02W = (-1) 2 Ra- — — — ■ — — — — ■ e 2 = 0,l, 5.10 

1 — a tanh(a; + 71) tanha(x + 72) 

(j-tanh(a; + 7i)tanha(a; + 72) 

03W = (-1) 3 ^ 7 — T7 — : — v— ( — : r, £3 = 0,1. 

1 — crtann(2; + 71) tanno"(j; + 72) 

5.3 The structure of the moduli space of kinks 

To describe the structure of the two-dimensional moduli space of solitary waves solutions (5.10) 
it is convenient to rely on a re-shuffling of the coordinates of the moduli space: 7 = —71, 
7 = 72 — 71 . 7 determines the "center" of the energy density of a given kink and 7 distinguishes 
between different kinks in the moduli space by selecting the kink orbit in S 2 , see Figs. 5 and 6. In 
the Figs. 3 and 4, however, we have plotted the components of the kink profiles for two different 
choices of 7 and 7 = 0. We remark that none of the three components are kink-shaped, the <j)\ 
and 4>3 components are bell-shaped and the </>2-components have a maximum and a minimum. 
The difference is that the x —> —x reflection symmetry is lost for 7 = 5. The profiles of 
these two solitary waves tend to the South pole in both x —> ±00. Therefore, these kinks are 
non-topological kinks living in Css- 

(7, 7) are good coordinates in the solitary wave moduli space characterized by the orbits 
in P2- The inverse mapping to S 2 is classified by the different choices of ex, €2 and €3: 

• The value of €3 selects the topological sector of the solution. Clearly: 

£3 = 0=^ lim 03 = —R, £3 = 1=^ lim ^3 = R 

x- > roc x— >±oo 

and the equations (5.10) determine two families of non-topological (NTK) kinks, one 
family living in Css, £3 = 0, and the other one belonging to Cnn, e 3 = 1- They are 
homoclinic trajectories of the analogous mechanical system. 
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Figure 3. Graphics of (5.10) corresponding to the following values of parameters: R = 1, er = 0.7, 
7 = 0, 7 = 0, i.e. the maximally symmetric generic kink (it is chosen e\ = £2 = £3 = 0). 




Figure 4. Graphics of (5.10) corresponding to the following values of parameters: R = 1, a = 0.7, 
7 = 0, 7 = 5, ex = e 2 = £3 = 0. 




Figure 5. NTK kink-orbits, corresponding to: R = 1, a = 0.7, £1 = £2 = £3 = and 7 = (left), 7 = 5 
(middle). Several kink orbits in C$s corresponding to several different values of 7 (right). 



• e\ = provides the NTK kink orbits running in the 4>\ > hemisphere and ei = 1 
corresponds to the NTK orbits passing through the other (</>i < 0) hemisphere of §> 2 . 

• £2, however, does not affect to the kink orbit and merely specifies the kink/antikink cha- 
racter of the solitary wave, i.e. the sense in which the orbit is traveled. 

In Fig. 5 (right) several NTK orbits are shown in §> 2 . The sector Css and the 4>\ > 
hemisphere have been selected. All of them cross the = meridian at the same point: 
F = (Ra, 0, Ra) - in sphero-conical coordinates this point is the low right corner Ai = A2 = o 1 
of P 2 ). Therefore the point F is a conjugate point 1 of the South pole S of S 2 . 

In Fig. 6 the graphics of the kink energy densities of two members of the NTK family are 
plotted. We remark that the energy density appears to be composed by two basic elements. 
This behavior is more evident for I7I 3> because in such a range the energy density resembles 
the superposition of the energy densities of the two embedded sine-Gordon kinks (compare with 
Fig. 2 (right)). All the NTK kinks, however, have the same energy that can be computed from 

1 We shall show in Section 7 that there exists a Jacobi field orthogonal to each NTK orbit that becomes zero 
at S and F, confirming thus that they are conjugate points. We shall also profit of this fact to show the instability 
of the NTK kinks by means of the Morse index theorem. 
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-10 
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Figure 6. Kink energy density for the NTK determined by the constants: R = 1, a = 0.7, and 7 = 7 = 
(left), 7 = and 7 = 5 (right). 

the formula (4.6): 

E NTK = 2uR 2 (l + a) = E Kl + E K . 2 . 

There is a kink energy sum rule: the energy of any NTK kink is exactly the the sum of the 
two sine-Gordon, K\ and K2, topological kink energies. These sum rules arises in every field 
theoretical model with mechanical analogue system which is Hamilton-Jacobi separable and has 
several unstable equilibrium points, see, e.g., [16]. Moreover, looking at Fig. 4 as a posterior 
picture to Fig. 3 in a movie that evolves with increasing 7, it is clear that we find almost a K\ 
anti-kink in the i>0 region and a K% kink in the iCO region when 7 — > 00. This means that 
the solitary waves (5.10) tend to a K\ kink/antikink plus a K2 antikink/kink when 7 — > ±00. 
Thus the combinations K\ + K2 (kink/antikink and antikink/kink) form the closure of the NTK 
moduli space and belong to the boundary. Alternatively, one can tell that there exist only two 
basic (topological) kink solutions and the rest of the kinks, the NTK family, are (nonlinear) 
combinations of the two basic solitary waves. 

6 The moduli space of § 3 -solitary waves 
6.1 Embedded N = 2 kinks 

In the massive nonlinear § 3 -sigma model we first count the three embedded sine-Gordon kinks: 
• The K\jK\ kinks and their anti- kinks living in the meridian 



{0? + 4 = r 2 } = {{h = h = 0} p| {4 +4+4+4=^ 



}} 



the intersection of the 4>2 = (f>s = plane with the S 3 sphere. 
• Idem for the K2/K2 and their anti- kinks: 



+ ^1 = r 2 } = {{^! = c/> 3 = 0} p| {4 + <& + 4>i + <& = r 2 



}} 



• Idem for the K3/K3 kinks and their anti-kinks: 



{4 +( f>i=R 2 }^ {{^ = <p 2 = 0} p| {4 + 0! +<%+ 4 = r 2 



}} 



The novelty is that all the non-topological kinks of the massive nonlinear S 2 -sigma model are 
embedded in the S 3 version in three different S 2 "meridian" sub-manifolds. 



Solitary Waves in Massive Nonlinear § -Sigma Models 



13 



1. Consider the two-dimensional S| sphere: 

{4 + 4 + 4 = R 2 } = {{0i = o} p {0? + 4 + 4 + 4 = R 2 }} , 

the intersection of the cj)\ = 3-hyperplane in M 4 with the § 3 -sphere. 

The restriction of the N = 3-model to these two-manifold collects all the solitary waves 
of the N = 2-model. We shall call NTKj kinks to the N = 2 NTK kinks of the N = 2 
model that lives in Sf. The K 2 /K% and K 3 jK\ sine-Gordon topological kinks/anti- kinks 
also live in in S 2 and belong to the boundary of the NTKj moduli space. 

2. Idem for the S 2 j sphere: 

{4 + 4 + 4 = R 2 } = {{0 2 = o} p| {4 + 4 + 4 + 4 = r 2 } } . 

The non-topological kinks running on this 2-sphere will be termed NTKjj. The sine- 
Gordon kinks/antikinks are: K\jK\ and K 3 jK% i . 

3. Idem for the S 2 n sphere: 

{4 + 4 + 4 = R 2 } = {{03 = o} n {4 + 4 + 4 + 4 = R 2 }}- 

The non-topological kinks running on this 2-sphere will be termed NT Km. The sine- 
Gordon kinks/antikinks are: K\/K\ and K 2 /K^. 

In the "equatorial" sphere 4 + 4 + 4 = R 2 > 04 = 0> °f however, there are no kinks because 
the equilibrium points, the North and South poles, are not included. 



6.2 Generic solitary waves 

Away from the meridian 2-spheres S|, S'jj and Sjjj considered in the previous sub-section, there 
is a three-parametric family of generic kinks. In N = 3 dimensional sphero-conical coordi- 
nates (4.2): 



k 2 



UlU 2 U 3 2 _ (cr| - Ul)(4 ~ U 2 ){(T2 - U 3 ) 



1 — u _2-2 ' 92 — u -1 



2) 

(gj ~ U l)(^3 ~ ~ Ug) ,2 (\ - U X ){\ ~ U 2 ){\ ~ U Z ) 
3—^0 2-21-2 =2~\ ' 04 _ U 2~~2 ' 



,2 



such that: < u\ < a\ < u 2 < erf < u 3 < 1, the § 3 -sphere of radius R is characterized by 
the equation: no = R 2 ■ The change of coordinates (6.1) corresponds to a sixteen-to-one map 
from the points of S 3 out of the coordinate 3-hyperplanes (j) a = 0, a = 1, . . . , 4, to the interior of 
the parallelepiped P3, in the u 2 , «3)-space, defined by the inequalities above. The meridian 
2-spheres are eight-to-one mapped into the boundary of P3, where some inequalities become 
identities: u\ = for §|, {u\ = g 2 }\J{u 2 = for §//, and {u 2 = a 2 } \J{u 3 = ^3} i n the 
S|jj case. The one-dimensional coordinate meridians of these 2-spheres, in particular those that 
are the embedded sine-Gordon kink orbits are four-to-one mapped into the edges of P3, see 
Fig. 7. Finally, the North and South poles of S 3 are two-to-one mapped into one vertex of P3: 
v ± = (ui,u 2 ,u 3 ) = (0, erf, erf). 

The quadratic polynomial V((f>i, (j) 2 , 4> 3 , $4) becomes linear in the sphero-conical coordinates: 



v 2 R 2 

V(ui, u 2 , u 3 ) = —— (ui + u 2 + u 3 -a 2 - 4) 
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Figure 7. Parallelepiped P3. The K\, K2 and K% sine-Gordon kink orbits are depicted (in blue) at the 
corresponding edges of P 3 . 



but we write it in "separable" form: 

• 2 R 2 f Ml (ui — of )(ui - of) 
(«i - u 2 )(u 1 - u 3 ) 



V(ui,U2,U 3 ) 



u 2 (u 2 - Q-j){u 2 - <j£) u 3 (u 3 - a 2 ){u 3 - o|) 
(u 2 - ui)(u 2 - u 3 ) (u 3 - Ui)(li 3 - u 2 ) 



The reason is that in this form it is obvious that the PDE (4.7) reduces to the system of three 
uncoupled ODE that can be immediately integrated to find (4.10): 

W( Ul ,u 2 ,u 3 ) = R 2 ((-l) 5 i7T^ + (-1)^71^ 

+ (-lY^VT^s) , 6 1 ,6 2 ,5 3 = 0,1. 



The Bogomolnyi splitting of the energy (4.6) prompts the first-order equations (4.11): 

dx (-l) Sa du a 
TTT? — \ = n — / ^2T7 -2\ i~\ a = 1,2, 3. (6.2) 

U'[U a ) 2u a (U a - a%)(u a - V 3 )\/l - U a 

It is convenient to pass to the variables: 

si = (-l)*vT=Hr, s 2 = (-l) 5 *VT^2~, s 3 = (-l) 53 vT^3~ 
because the quadratures in (6.2) become: 

Eds a , \ v ds a \ -» ds a . . 
5- = -dx, } —5 5 = -dx, } —5 77 = -dx, (6.3) 

a=l " a=l z a a=l J a 

where only enter rational functions. The integration of (6.3) is straightforward: 

arctanh si + arctanh s 2 + arctanh S3 = — (x + 71), 

-o 2 (x + 7 2 ), (6.4) 

-cr 3 (x + 7 3 ). 

The addition formulas for the inverse hyperbolic functions, and the definition of "Vieta variab- 
les": 



arccoth 


si 


+ arctanh 


S2 


+ arctanh 


S3 




(72 




09 




cr 2 


arccoth 


si 


+ arccoth 


S2 


+ arctanh 


S3 




03 




03 




03 



^4 = si + s 2 + s 3 , -B = sis 2 + S1S3 + s 2 s 3 , C = S1S2S3 
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reduce the system of transcendent equations (6.4) to the linear system: 

A - tiB + C = ti, 
a\t 2 A - a 2 B + t 2 C = el, 
a\A - a 3 t 3 B + C = oJt 3 , 



where: t\ = tanh( — {x + 71)), t 2 = tanh(— a 2 (x + 72)), *3 = tanh(— a%{x + 73)). From the 
solutions A, B and C of this linear system we obtain the s±, s 2 and S3 variables as the roots of 
the cubic equation: 



„3 A„2 



As z + Bs-C = 0. 



The standard Cardano parameters 



B A 2 1 A 3 

q = — — , r = -(3C — AB) + — — , 6 = arccos ■ 



3 9' 6 V ' 27' y^gS 

provide the three roots si, s 2 and S3 by Cardano's formulas, and recalling that: u a = 1 — s 2 ., 
a = 1, 2, 3 we finally obtain: 

A 



Ul (x) = 1 - ( — + 2 v /= gcos - 

n 2 i •'• ) - 1 - ( ^ + V^q (- c °s I - \/3 sin I ) ) . (6.5 ) 

l'A , / 9 r 9 

"3 1 ' ) - 1 - ( — + V-q I - cos - + V3 sm - 

which in turn can be mapped back to cartesian coordinates. The explicit analytical expressions 
for the generic kinks are: 

, , Tseu,- - (°"2 -0-3*2*3) sech(x + 71) 

^l(x) = (-1) 1 RtT2<T3^2 2777 =27^' 

cr 2 o-| - (erf - og)*i* 2 - o 3 crf * 2 *3 
1 f \ , 1 \€o t~>— l~2 2 (*1 ~ 0-3*3) sech<r 2 (x + 72) 

<h{x) = (-l) e2 i?0 2 A/al - O^ — 7"2 2777 T2— ' 

O2O3 - (Cr 2 - Cr 3 )M 2 - 030-2*2*3 , . 

1 f \ 1 1WD - ^1 2" (°2 -*i*2)sechcr 3 (x + 7 3 ) 



bAx) = (-IY*R 



a 2 a\ - (of - cr|)*i* 2 - o 3 o\t 2 t 3 ' 
-ct|ct 3 *i - (of - of - 0-203*1*2) * 3 
cr 2 of - (cr| - cj|)*i*2 - o 3 of *2*3 



where e a = 0, 1, V a = 1, . . . , 4. Equations (6.6) define a three-parameter family of kinks in the 
S 3 -sphere. 



6.3 The structure of the moduli space of kinks 

In the Figs. 8 and 9 two generic kink profiles have been plotted for representative values of the 
parameters. 

A glance to the plots of 04 as function of x shows that the kinks are topological interpolating 
between the North and the South poles. In fact, all the generic kinks are topological. The ana- 
lysis of the asymptotic behavior of the fourth field component in (6.6) determines the topological 
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character of all the solitary waves solutions. The generic kinks belong to either the Cns or the 
Csn sectors of the configuration space. The limits 

lim (Mx), fa(x),<h(*)> Mrf) = (0, 0, 0, (-l) e4 i?), 

x— >— oo 

lim (Mx),Mx),Mx),Mx)) = (0,0,0,-(-l) e4 i?) 

x— Yao 

show that €4 = implies that the TK kink family lies in CnSi whereas = 1 sends the 
TK family to live in Csn- The other sign options, (— 1) €1 , (— l) e2 and (— 1) £3 , determine the 
three-dimensional hemispheres accommodating the different TK kink orbits, as well as the 
kink/antikink character of the solutions, in a smooth generalization of the N = 2 NTK kinks. 

It is possible to draw the kink orbits in the P3 parallelepiped using the solutions (6.5) in 
sphero-conical coordinates. In the Fig. 10 (right) several TK orbits are drawn together. The 
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Figure 10. TK orbit corresponding to: 71 = 72 = 73 = (left). The case: 71 = 0.8, 72 = 2 and 73 = 2 
(middle). Three orbits corresponding to different values of the parameters (right). 




Figure 11. Kink energy densities for different values of the constants 71, 72 and 73. 



energy of all the TK kinks is: 

E TK = 2uR 2 (l + a 2 + <7 3 ) = E Kl + E K2 + E Kz . 

There is a new kink energy sum rule specifically arising in the N = 3 model. In the Fig. 11 
it can be viewed how the energy density reflects this energy sum rule: the energy density of a 
given member of the TK family is a "composition" of the energy densities of the sine-Gordon 
kinks K\, K2 and if 3, see Fig. 11 (middle) and (right). This structure is completely similar 
to the kink manifold of the field theoretical model with target space flat M 3 and the Gamier 
system as mechanical analogue problem, see [12] and [16]. In fact, the equations are identical, 
even thought in that case the problem is separable using elliptical coordinates in M 3 , [12, 16]. 
The calculation of the limits of the TK family when the parameters 71 , 72 and 73 tend to ±00 
is not a simple task, but it is possible to show that all the possibilities arising from the energy 
sum rule are obtained for the different limits (see [16] for details; the model is different but the 
equations are the same and thus the analytical calculation of the limiting cases coincides). In the 
Fig. 10 (right), for instance, a TK kink orbit (in red color) close to the limiting combination if 3 
plus a member of the NTKjjj family, is plotted. 

Finally, the structure of the kink manifolds for the N = 2 and N = 3 models allows us 
to foresee that in higher N, the generic kinks will be topological if N is an odd number, and 
non-topological for even N. In fact, the kink energy sum rule will establish that the energy of 
a generic solitary wave is the sum of the energies of the N sine-Gordon kinks embedded in the 
model. 

7 Kink stability 

The analysis of the kink stability requires the study of the small fluctuations around kinks. For 
simplicity, we will only present the explicit results for the topological kinks of the § 2 model, and 
also a geometrical analysis in terms of Jacobi fields that allows us to demonstrate the instability 
of the non-topological kinks. The generalization to the arbitrary N case is straightforward. 
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7.1 Small fluctuations on topological kinks, K 1 and K 2 

Taking into account the explicit expressions of the K\ and K2 sine-Gordon kinks of the § 2 model, 
we will use spherical coordinates in S 2 in order to alleviate the notation. Thus we introduce the 
coordinates: 



iisin^costf 2 , 



Rsin6 1 sm9 2 , (f> 3 = Rcos9 1 



with 6 1 6 [0,7r], 9 2 £ [0,27r). The metric tensor over S 2 is written now as: ds 2 = R 2 d6 1 d6 1 + 
R 2 sin 2 6 l d6 2 d6 2 . The associated Christoffel symbols and Riemann curvature tensor components 
will be: 



1 



L 22 - 

-^212 



sin20\ 



12 



T21 = cotan^ 1 , 



"^122 



sin 



2 Q 1 



p2 
^121 



~ n 211 



The analysis of small fluctuations around a kink 9k{x) = (6] < (x), 2 K (x)) is determined by 
the second-order Hessian operator: 



AkV = ~ (V^ V^r? + R(9' K , r,)9' K + V„grad V) , 



(7.1) 



i.e. the geodesic deviation operator plus the Hessian of the potential energy density, valuated at 
the kink. r\ denotes the perturbation around the kink. Let 9{x) denote the deformed trajectory, 
9(x) = 6k{x) +t](x), with r](x) = (rj 1 (x) , r] 2 (x)) , then we introduce the following contra-variant 
vector fields along the kink trajectory, 77, d' K G T(TS 2 \k )- 



r ? (x) = r/ 1 (x)^ T +r ? 2 (x)^2, 



9 k W + 0k 



d 



We will use standard notation for covariant derivatives and Riemann tensor: 

V^r? = ( V '\x) + T) k r?0'k)^ R(9' K , V )9' K = 9%r? (x)^ k ^ 
The geodesic deviation operator and the Hessian of the potential read: 

; ,; + r(6' k , v )e' K = v* v e >n + RK, v)o' K , 



dx 2 

Vt] grad V 



V 



d 2 v 



dV 



nil. 



T ij QQk ) 9 J QQl 



evaluated at 9k{x). In sum, second-order kink fluctuations are determined by the operator: 



A K T] 



^- + R(e' K , V )9' K + V v gTadV 

d 2 n l r 1 

- cos 26 l K y (6%) + a 2 + a 2 cos 2 8 2 K ^ r/ 1 - sin 29^6% 

(l + cos2^)« + ^ 

2cotan#i<9£^- + ( cotan - d'hB'lW + 

dx 



j//2 -2 



sin 29 2 K 



V 

d 2 V 2 
dx 2 



drj 2 
dx 

d 
W 



+ 2cotan^^ + (cotan^ 1 - (tf£) 2 - cos 2 ^(^) 2 ) N ° 



de 2 ' 



(7.2) 
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The spectrum of small fluctuations around K2/K2 kinks. The K2/K2 kink solu- 
tions (5.1), (5.2) are written, in spherical coordinates, as follows: 



4 2 = 2arctane ±CT ( :c - a;o \ 



2arctane ±ff(a: - Xo) , 



°K 2 



9 2 



7T 

2 : 
3vr 

T' 



where the (±) sign determine the kink/antikink choice, and we will take xq = for simplicity. 
Plugging these K2/K2 kink solutions in (7.2), we obtain the differential operator acting on the 
second-order fluctuation operator around the KijK\ kinks: 



dx 2 



2a 2 



cosh ax 



_d_ 

de 1 



+ 



d V , drr _ 2 2 

+ 2atanhcjx— h a r\ 

dx z dx 



d 



.(7.3) 



We know from classical differential geometry that this expression can be simplified if one uses 
a parallel basis along the kink trajectory. Therefore let us consider the equations of parallel 
transport along K2 kinks: v = 0, where v is the vector-field: v(x) = v^(x)J^ + v 2 {x)-^ I . 
These equations are written explicitly as: 



dv l 
dx 



' dy^_ 
dx 
dv 2 







cotan(2 arctan e ax ) 9 

, + a ^-j V = 

ax cosn ax 



Therefore the vector-fields: v\ = -J^ and V2(x) = coshax-^p form a frame {^1,^2} in Y(T% 2 \k 2 ) 
parallel to the K2 kink in which (7.3) reads as follows: 



v (x) 
v 2 (x) 



cosh ax. 



^K 2 V 



dx 2 



2a 2 



cosh 2 ax 



>1 



v\ + 



'2^2 



d 2 7] 

dx 2 



+ 1 



2a 2 



cosh 2 ax 



V 



V2, 



where r\ = fj vi+fj 2 V2, t] 1 = 77 1 , and rj 2 = cosh axff . Thus the second-order fluctuation operator, 
expressed in the parallel frame, is a diagonal matrix of transparent Posch-Teller Schrodinger 
operators with very well known spectra. In fact, in the v\ = ^ direction we find the Schrodinger 
operator governing sine-Gordon kink fluctuations, as could be expected. The novelty is that we 
find another Posch-Teller potential of the same type in the vi = gt^ direction, orthogonal to 
the kink trajectory. The spectra is given by: 

• v\ direction. There is a bound state of zero eigenvalue and a continuous family of positive 
eigenfunctions: 



t)q{x) = sech ax, sjp = 0, 

fjl( x ) = e ikax (t&nhax - ik), e (1) (fc) = a 2 (k 2 + l) . 



coshers-^ direction. The spectrum is similar but the bound state corresponds to 



V2 

a positive eigenvalue: 



section, e 



(2) 



1 - a z > 0, 



fg( x ) = e ikax (t<mhax - ik), e i2) (k) = a 2 k 2 + 1. 
Because there are no fluctuations of negative eigenvalue, the K2/K2 kinks are stable. 
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The spectrum of small fluctuations around K\/K^ kinks. A similar procedure shows 
that the K\ kink/antikink are unstable. The K\/K* kink solutions (5.1), (5.2), in spherical 
coordinates, are as follows: 



2arctane ±(:r - :Co) , 
\» =2arctane ±(a; - a;o) , 



9 2 



0. 

: 7T, 



where we find again that the (±) sign determine the kink/antikink choice, and it will be taken 
xo = for simplicity. By inserting the K\jK\ solutions in (7.1) the second-order fluctuation 
operator around the K\/K* kinks is found: 



A Kl r] = A K *t] 



dx 2 



+ 1 



cosh x 



d 

w 



drr _ 2 2 

— y + 2 tanhx— a r\ 

dx z dx 



d 

w 



(7.4) 



An analogous calculation lead to the parallel frame along the K\/K* kink orbits: 

d 



{ Wl ,w 2 }€r(m 2 \ Ki ), 

And thus (7.4) becomes: 



u>i 



d 



W2(x) = coshx 



de 2 ' 



A Kl 7] = A K *rj 



+ 



dx 2 
_d 2 rf 
dx 2 



+ 1 



cosh x 



v 1 



+ a 



cosh 2 x 



f] 



w 2 



(7.5) 



fj , T] 2 = cosh xfj 2 . 



with 77 = ^wi + rj 2 W2, T} 1 

Again, the second-order fluctuation operator (7.5) is a diagonal matrix of transparent Posch- 
Teller operators. In this case, there is a bound state of zero eigenvalue and a continuous family 
of positive eigenfunctions starting at the threshold £^(0) = 1 in the w\ = g|x direction: 



f)l(x) 



1 



.(1) 



cosh x ' 



fjl( x ) = e ifca; (tanhx - ik), e^(k) = (k 2 + l) . 



In the W2(x) = cosh x JL direction, the spectrum is similar but the eigenvalue of the bound state 



is negative, whereas the threshold of this branch of the continuous spectrum is £^(0) = a 2 : 



€2-1 0) 



1 

cosh x ' 



a 2 - 1< 0, 



(tanhx — ik), e^ 2 \k) 



k 2 + a 2 . 



Therefore, Ki/K* kinks are unstable and a Jacobi field for k = ia arises: fjj(x) = e CT:r (tanhx — 

Finally, it is remarkable that this results can be reproduced without difficulty in the cases 
> 2 by using hyper-spherical coordinates in S . The only stable kink that is present in the 
model is the sine-Gordon Kjy kink, that corresponds to the minimum of the kink energies. 



7.2 The stability of the NTK kinks 

The analysis of the spectrum of small fluctuations around the kinks of the NTK family is not 
an easy task, and thus we will proceed by another way, see [16, 19]. Having obtained explicit 
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Figure 12. Components of the Jacobi fields (7.6) for a — 0.7, ei = 62 = £3 = 0, and 7 = 7 = (left), 
7 = 1,7 = 2 (right). 



expressions § NTK (x; 7, 7) for the NTK family of solutions (5.10), we know from classical Morse- 
Jacobi theory about conjugate points, that the vector-fields ^Wr~ an d d ^g^ K are Jacobi fields 
along the NTK trajectories, i.e. zero-modes of the spectrum of small fluctuations around the 
solutions. In fact, taking into account that the 7 parameter simply determines the "center" of 
the kink, it is easy to conclude that ^ is tangent to the kink trajectories, and thus only 9 ®q^ K 
is a genuine (i.e. orthogonal to the kink-orbits) Jacobi field for this family of solutions [16]. 
By derivation on (5.10) with respect to 7 parameter, it is obtained the Jacobi field: 

tNTK -\ d& Raa sech(cr(a; — 7 + 7)) tanh(x — 7) 

J (x; 7, 7) =- - 



<9 7 (1 + a tanh(x — 7) tanh(c(x — 7 + j)) 2 

( -9 
x (— l) ei cj sech(x — 7) sech(cr(x — 7 + 7))^— — 

d 

+ (-1)' 2 (cr tanh(x - 7) + tanh(cr(x - 7 + 7))) — — 



d 

(-l) £3 asech(CT(x-7 + 7))— ) . (7.6) 
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By a direct calculation, it can be checked that 

lim J NTK (x; 1 ,j) = Q, J NTK ( T ,l,l) = 

x— >±oo 

and thus the point x = 7 in the kink trajectory, i.e. $ (757, 7) = ((— l) €l Ro; 0, (— l) e:1 Ra) 
is a conjugate point of the corresponding minima S (or N depending on the choice of 63). In 
Fig. 12 there are depicted the components of (7.6) for several different cases. 

The existence of a conjugate point establishes, according with Morse theory, that these kink 
solutions are not stable under small perturbations. In fact, the Morse index theorem states 
that the number of negative eigenvalues of the second order fluctuation operator around a given 
orbit is equal to the number of conjugate points crossed by the orbit [19]. The reason is that 
the spectrum of the Schrodinger operator has in this case an eigenfunction with as many nodes 
as the Morse index, the Jacobi field, whereas the ground state has no nodes. The Jacobi fields 
of the NTK orbits cross one conjugate point, their Morse index is one, and the NTK kinks are 
unstable. 

Finally, it is possible to extend these results to the N = 3 case, where two Jacobi fields 
appear, and the instability of the generic TK family is showed. The procedure is the same, but 
the extension of the expressions is considerably bigger, thus we will not include this calculation 
here. 
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